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Abstract 

We consider an open bipartite quantum system with dissipative dynamics generated 
by Lie = Lo+eLi, where Lo,i are generators of Lindblad type and < e « 1. In order 
to study the entanglement of the stationary states of L^, we develop a perturbative 
approach and apply it to the physically significant case when Lq generates a reversible 
unitary dynamics, while Li is a purely dissipative perturbation. 

1 Introduction 

An open quantum system dynamics of Lindblad type [U |2] is an effective description of the 
action of an environment E on a quantum system S weakly coupled to it. In general, the 
environment acts as a source of dissipation and noise; in spite of this, decoherence is not 
the only possible consequence. If suitably engineered, the coupling with the environment 
may also generate coherence and even entanglement [3], this possibility depending on the 
trade-off between dissipative effects and environment induced mixing [U El El El El El HH] - 

Of particular interest is under which conditions the presence of an environment may induce 
convergence to asymptotic states with definite entanglement properties [HI |121 HSl [H] . In 
fact, controlling the coupling to the environment could then be used for preparing states 
with definite entanglement content From this viewpoint, it is of great importance to 
know 1) the invariant states of a given Lindblad dynamics and 2) whether any initial state 
converges asymptotically to some stationary state. A part from some older [161 El [E] cind 
more recent results [191 EHl EU ES] , a full characterization of the asymptotic properties of 
open quantum systems and their asymptotic behavior is still to be achieved. 

In the following, we will focus upon the following scenario: consider two finite-level systems 
and 5*2, not directly interacting with each other, whose reversible, unitary dynamics is 
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generated by a Hamiltonian H = Hi + H2 via the generator Lo[p] = —i[H,p]. If weakly 
coupled to a same environment, on a long time-scale, they undergo an open, dissipative 
dynamics generated by = Lq + £ Li, where Li is a generator of Lindblad form and e 
measures the weakness of the coupling to the environment. In general, the addition of the 
perturbation term e Li diminishes the number of invariant states with respect to those of 
Lq; however, in finite dimension, at least one invariant state will always survive and, by 
continuity, will be close to them: the issue is whether the remaining invariant states may be 
entangled or not. 

Suppose the spectrum of the Hamiltonian H = H1+H2 be non-degenerate, then Lo[po] = 
only if p is a separable state. Intuitively, if such states are well inside the closed convex subset 
of separable states, no dissipative perturbation Li could provide entangled states such 
that lLe[pe] = 0. Indeed, by continuity, such asymptotic states are perturbations of those of 
Lo, namely p^ = po+e pi+o{e), and thus remain separable if po is separable. On the contrary, 
for separable stationary states po on the boundary of the subset of separable states, it should 
be possible to construct entangled ps by suitably engineered, small dissipative perturbations. 

In the following, we give mathematical ground to these expectations by developing a 
systematic perturbation expansion of the states p^ that are invariant under generators of the 
form Le = Lo + ehi where Lq and Li are generic Lindblad type generators. 

2 Perturbation theory 

Let 5* denote a li-level system with observables X from the full matrix algebra Md{C) and 
states (density matrices) p from the convex subset S{S) C Md(C) of positive matrices of 
trace 1. If is weakly coupled to its environment E, its time-evolution is conveniently 
approximated by a Markovian Lindblad-type dynamics: p ^ pt = Jtlp] = exp (tlL)[p] [H |2] 
where L is the generator of the semigroup of trace-preserving, completely positive maps 7^, 
7^07^ = 'jt+s- It incorporates in an effective manner the noise and dissipation due to the 
environment via a master equation of the form 

dtpit) = L[p(t)] = -t[H , Pit)] + D[p(t)] , (1) 

where Md{€) 3 H = W , while 

B[p] = Y,{h^phi - l{hih^,p}) , (2) 

a 

where , Ea ^iK e Md{€). 

We shall denote by 5^ C 5 the subset of stationary states of 7^: they satisfy L[p] = and 
form a convex subset of /C(L), the kernel of the generator. 

The time-evolution generated by affects the states of the system, while its observables 
evolve according to the semigroup of dual maps 7^"" : X ^-> X(t) = -fJ[X] generated by 

dtX{t) = LT[X(t)] =t[H, X{t)] + DT[p(t)] , (3) 
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where 

mX] = J2{hlXh^ - ^{hih^^x]) . (4) 

a 

About the structure of the 7t-invariant states, we have [16] 
Proposition 1. The time-average 

X^G^[X]= hm 1 rdtjJlX] (5) 

is a well-defined unital, completely positive map; its dual map G \ S ^ defined by 

Tr(G[p]X^ = Tr(^pG^[X]) VX e Md(C) , p G 5 (6) 

is a completely positive, trace preserving map which associates a state p e S to a 'jt-invariant 
state: LoG[p] =0. 

If jt possesses a faithful invariant state, that is a full-rank density matrix p* = 7t[p*], 
then the -invariant observables X = jJ[X] form a subalgebra C Mrf(C) and is a 
conditional expectation onto M^: 

G^[Yi X Y2] = Fi G^[X] Y2 V Fi,2 e , X G Md(C) . (7) 

ControUing the structure of and, in particular, whether an initial state p G 5 converges 
to G [p] is a complicated matter with some partial clues [lU [T71 |T8] . Recently, such an issue 
has become again subject of study [HI EOl EB 123] also because of its increasing importance in 
quantum information [15]. Concerning 7t-invariant states, the following result was obtained 
in [22] (see also [23])- 

Proposition 2. Let L be the generator of a Lindblad-type dynamics 7t; one can always 
construct orthogonal stationary states pj of 'jt: ^[Pj] = and pjp^ = unless j = k. 



The following ones are instances of some possible scenarios. 



Example 1. Let L[p] = ~i[H,p], where the Hamiltonian H = Ylj=i \ j) Ul ^ non- 
degenerate spectrum; then, the stationary states pj of Proposition [H are the orthogonal one- 
dimensional eigen-projectors \j){j\- For later application, we need extend the map G of 
PropositionUl to the whole matrix algebra M(i(C); in the present case it reads 

d 

G[X] = Y,{j\X\j)\j){j\ . (8) 

i=i 

Clearly, because of the oscillatory behavior, there is no tendency to equilibrium: 'jt[p] does 
not converge to G[p] when t -\-oo. 



3 



Example 2. Let H = in (QJj and ha = in where {\a)} is an orthonormal basis 

o/C"'. Then, 



^p] = \^P){^|J\ - p^^,[p]=e-'p+(l-e-')\^){^P\ . (9) 

Hence, p = |'?/')('?/'| is the only stationary state and all others converge to it asymptotically: 
7t[p] I — > The corresponding map (G is given by (on Mrf(C)j 

G[X] = Tr(X) |^)(^| . (10) 

Example 3. Let if = m (QP and hi = \il>){l\ in where WipW = 1, {|«)}a=i is an 
orthonormal basis and < \ipi\ = |(1|V')| < 1- Then, 

I^[P] = (1|P|1)I^)(^|-^|1)(1|P-^P|1)(1| • (11) 

Setting X^p = {a\X\(3) for all X e Md{€), p = 1 - {ipil^ > 0, ^„ = {a\ip), a ^ 1, the 
equations, 

Pll = -PPII , Pla = V^lCPll - 2 ' ^ i'a'ip^Pu 

can easily be solved yielding 

Pu{t) = Pll , Pia(t) = e"*/^ + V'iCPii ^ _ 
1 — e^'^* 

Pa^{t) = Pa(3 H Plllpa^} , 

P 

where a, b in the second expression are fixed by the initial conditions. Then, 

G[X] = J2 + —^ar^i) • (12) 

a,l3>2 ^ 

All states such that p = QpQ, where Q = J2a>2 ^'^^ ''yt-invariant; also, '~it[P\ ' — > *G[p] 



Despite the abstract characterizations of [2T|, [23] , the convex subset of stationary states is 
difficult to control in practice; we shall thus concentrate on understanding how the invariant 
states of a semigroup '-j'f'' are modified by a perturbation Li of its Lindblad generator Lq. 
Concretely, we will investigate the set 5^ of stationary states of Lindblad-type dynamics 7^^'' 
generated by = Lq + £lLi, < e <C 1. By switching on the perturbation, the dimension 
of Se decreases, but. Proposition [1] ensures the existence of at least one stationary state. 

Lemma 1. Consider the generator = Lq + s'Li, where both Lo,i are generators in 
Lindblad form and the semigroups 7^°^ and 7^'^'' generated by Lq and = Lq + e Li . Let n(0) 
be the number of "^f^ -invariant orthogonal density matrices and n{e) that of '-yf^ -invariant 
orthogonal density matrices for < £ <C 1; then n(0) > n{e). 
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Proof: From Proposition O one can always choose density matrices such that Le[pj(£:)] = 
and TT{pj{e)pk{e)) = for j ^ k. In finite dimension, eigenvalues and eigen-projectors are 
continuous in e; therefore, should Pj{e) ^ Pki^) merge as e — )■ 0, the continuity of the 
Hilbert-Schmidt scalar product would be violated. 

Because of finite dimensionality, the solutions can always be expressed as converging series 
in powers of e 

p, = ^£"p„, (13) 

n>0 

where the operators pn must solve the iterative procedure 

(Lo + eLi) [p(£)] = Lo[po] + 5^£"(Lo[p„] + Li[p„_i]) =0 , whence (14) 

n=l 

Lo[po] = , Lo[p„] = -Li[p„_i] n > 1 , (15) 

where po is a stationary state of 7^^°^ Also, since Tr(p(£)) = 1, it follows that Tr(p„) must 
vanish at all orders. In the following, we discuss when p„ = — Lq ^[Li[p„_i]] are acceptable 
solutions. 

Definition 1. Let F = id — (G, where G is as in PropositionU^ since (G is trace-preserving, 
the image of Md{(C) by F consists of traceless matrices: Tr(F[X]) = for all X G M(i(C). 



Lemma 2. L ^ can he defined as a map from F[Mrf(C)] into itself. 

Proof: Notice that G, as a time-average, maps into the kernel of L and leaves it invariant; 
thus, from X = G[X] + F[X], it follows that L-^[X] is well defined on Md{€) 3X^0 
only if G[X] = 0. Then, is constructed as a linear map from the range of F into itself 
such that L o L"^ = o L = id on F(Mrf(C)). This guarantees that L"i[0] = 0; indeed, 
consider Z = h-^[X] - L-^[Y], with X = h[V] = Y = L[V + W],W ^ 0, h[W] = 0; then, 
F[Z] = F o L-i[X] - F o L~i[r] = 0. 

Example 4. In the case of Example Ui where L[p] = —i[H, p] and H is non-degenerate, 
G[X] = if and only if {j\X\j) = for all j ; one thus gets 

^-^M = ^E#^b)(^i- (16) 

Example 5. In the case of Example\^ ^[^] = ^ if o^^'^d only if Ti{X) = 0; one can verify 
that, on traceless matrices, 

h~^[X] = -X. (17) 
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Example 6. Finally, in Example\^ = and only if X is of the form 

X = Xn|l)(l| + 5^(Xi,|l)(a|+X,i|a)(l|) -^Q ^ V'.V'; 

where the only free entries are Xi^ and X^i, a > 1. Then, 
L-i[X] = --(Xn|l)(l|+ 5^ 

a>2 \ ^ ^ 

From the previous Lemma, it follows that, in order to solve for p„ in fllSp by inverting 
Lq, one has to ensure that (Go[lLi[pn-i]] = for all n > 2. The following Lemma gives a 
sufficient condition for this to be true. 

Lemma 3. Given = Lq + s\ji, z/Lo[po] = for a unique state po O'nd Li[po] 7^ 0, then 
Le[pe] = for a unique pe given by 

oo ^ 
n=0 u 



Proof: As Go maps into the kernel of Lq and is trace preserving, from the hypothesis of 

the lemma it follows that Go o Li[po] = A po- Then, A = Tr^Go o Li[po] j = Tr ^Li[po] j = 

implies Go o Li[po] = so that Lo can be inverted on Li[po] and one can solve the first 
recursive relation in (fTB]) . As Lg ^ maps into Fo[Mrf(C)] where Fq = 1 — Go and Gq is the 
trace-preserving map in corresponding to Lo, then Tr(pi) = 0. Iterating this argument 
yields the result. 

Example 7. For Lq as in Example [H there is only one invariant state so that Lemma 
applies. Furthermore, since Ijq'^[X] = —X / UP)) yields pe = ^1 — eLi j [p]. In such a case 

of a unique invariant state under Lo, we can make some preliminary considerations about 
the entanglement of the unique state invariant under Lo + eLi. Consider a separable pure 
state p = P ® Q ^ Mrf(C), where P = and Q = \x){x\; then, suitable non-local 

perturbations, Li may entangle it. Indeed, by partial transposition J^, p^ ^ p^ , operated on 
the second party with respect to an orthonormal basis starting with \x) , one gets 



P^Q + e(hi[P^Q]y + o{e) . (20) 
By projecting with Il± onto a subspace orthogonal to P ® Q, it follows that 

Tr(pf n^) =£Tr(n^(Li[P®Q])^) + o{e) . 



//Li[p] = —i[Hi ® 1 + 1 ® H2 + Hi2,p], where H12 is a non-local coupling of the two sub- 
systems, then the quantity 

Tr (p, n^) ~ -le Tr (u^([Hu , P ® Q]) ^) 

can be made negative by suitably choosing H12; then, one violates the positivity of partial 
transposition at order e and pe is entangled at that order. 

Entanglement can also be obtained via a purely dissipative time- evolution as the one gen- 
erated by Li as in /[TJ]) : indeed, choosing |1)(1| = P ®Q yields 

Tr(pfnx):^£Tr(nx(|^)(^|)^) , 

which can become negative by a suitable choice of entangled and Il±. 

The possibility of generating entanglement in the above two cases comes from the fact 
that the 0-th order state Pi ® P2 is on the border of the closed subset of separable states 
and can thus be moved into the open complementary subset of entangled states by suitable 
terms of order e. 

2.1 dim(ker(Lo)) > 2 

If, as in Examples [1] and |3l the kernel of Lq contains more than one stationary state, still 
one may seek a po such that Lo[po] = and 

(Go o ILi [po] = £1 [po] = , where Li := (Go o Li o Go , (21) 

so that the first order correction can be obtained as 

pi = -Lq^ o Li[po] . 

In order to continue the iteration in ([T^ and get 

P2 = -Lq^ o Li[pi] , 

again by inverting Lo, one has first to ensure that 

Go o Li [pi] = -Go o Li o Lo ' o Li [po] = , (22) 

and, analogously, for the higher order contributions to (!T3|) . 

Example 8. Consider the case where Lo,i[p] = — z[i/o,i 1 p\ with Hq non-degenerate. With 
Hq\E^-) = E^\E^) and using ^, one gets' 
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Hi, \E^){E^\ \E') \E^){E^\=0 



for all p. Then, with po = Ylt=oPk\^k){^k\! using IfT^) . one computes 

p, = -Lo ^ o U[po] = - E 1^ {E^mE',) \E^){El\ . 

From non- degenerate perturbation theory, the perturbation of \E^) to first order in e is the 
eigenvector of = Hq + e Hi given by 

Thus one sees that, to orders, reproduces pi with po = \E^){E^\. Furthermore, 

GooL,h]=.5^|^ \{ElmE^)\' {\El){El\ - |i?°)(E°|)=0. 

j^k J 

Thus, p2 = ^Lq ^ ° -''^i) [Po] Oj'i^d so on with higher orders. 
Unlike in the previous example, it may happen that ( 122|) is not satisfied by the chosen po- 

Example 9. Consider Lq as in ExampleUl and Li as in Example\^ the solution to 

d d 

LoM = and to U[p] = J] p,, (2|Li[|j)(j|]|z) = ^ - p,,) = 

i,j=l 1=1 

must have the form po = X]^=i 1^0) P Ij)(jT Then, a natural candidate for the first order 
perturbation contribution pi is, using Lq ^ as in Example\^ 

Pi = -i^o' ° U[po] = -Lo^ [|^)(^| - Pol =-^Yl ^^-^y^^^ \j){k\ . (23) 



Ea — El, 
j^k ^ 



However, with this choice, it turns out that 



Go o Li[pi] = Go [|^)(^| - PiJ = E ^ • 

A possible strategy to overcome the problem exposed by the previous example is as follows: 
in Lemma El Lq ^ is defined as a map from the range of Fq into itself. Thus, given a first 
order perturbation contribution pi = — Lq ^ o Li[po], one can always add to it (see the proof 
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of the Lemma) 0"i G Mrf(C) such that Lo[cri] = whence Lo[pi + cri] = Lo[pi] = — Li[po]- 
One can thus try to find an appropriate 7^'^'' -invariant matrix ai such that 

Go o 1^1 [Pi] + Go o Li[ai] = Go o Li[pi] + Li[ai] = , 

where we have used that Go[o"i] = cti. Thus, if such ai can be found it is of the form 

ai = -L^i o Go o Li [pi] = o Gq o Li o Lg ^ o Li [po] , (24) 

where the inverse L^^ of Li is defined as in Lemma[2]and thus Tr((Ti) = 0. Then, one would 
obtain the second order perturbation contribution p2 = — IL,^^ ° ^i[pi + 

Remark 1. Of course, the existence of ai is equivalent to the invertibility of hi. In gen- 
eral, Li is not a generator of Lindblad form; namely, it does not generate a semigroup of 
completely positive maps on the set of all matrices, even if Li does. However, in the next 
section, we shall consider the setting of Example and prove that Li generate a positivity 
preserving semigroup on the density matrices commuting with H . 

3 Dissipative Perturbation of a Unitary Evolution 

In the following, we restrict to a less general situation than the ones addressed in the previous 
section; namely, we will stick to purely dissipative perturbations Li = D as in ([2]) of a 
generator Lo [p] = ~i[H p] as in Example [1] 

L,[p] = -i[H,p\ + £D[p] , D[p] = [Kphi - \{hlh^,p]) . (25) 

a 

Lemma 4. The map D = Gq o D o Gq, with Go given by generates a positive, trace 
preserving map on the '-yf'' -invariant states p that commute with H . 

Proof: If 5o 9 p = E?=iPii 

d d 

fi[p] = E b)oi = E E {\(3K\^)\' - s., {j\hih^\j)) . (26) 

i,j=l «J=1 " 

Then, 

2 

Pn = ^PiiYl - {j\hlha\j)^ > -hpjj , (27) 



i=l 



where h 



• Therefore, the eigenvalues of any p G Sq remain positive while 



evolving with exp(tD). 



9 



Example 10. Consider Li as in Example\^ on density matrices Sq 3 p = J2j=iPjj 

that commute with H, one finds that exp(i:D)[p] converges to the unique invariant state 

Po = EU Indeed, ^ yields 

Pn = - Pn . pM = I V^(J) I' (l - e"*) + Pi^e"* . 

As already observed, even if one knows the structure of the invariant states of L(£:), it 
remains to be proved that one actually has asymptotic convergence to them. As showed 
in Lemma [H even a very weak perturbation eD of Lq in general decreases the dimension 
of the kernel of Lq; this is why adding a suitable engineered dissipative perturbation can 
be used to drive a system into a certain stationary state which may even be chosen to be 
pure However, one must check that all other eigenvalues of get a negative real 

part. This cannot hold in general [23]: purely imaginary eigenvalues can remain, but only 
if the Lindblad dynamics becomes trivial when reduced to the subspace which supports the 
stationary states. The following Lemma provides sufficient conditions for this not to be the 
case. 



Lemma 5. Assume that, given the generator ^ME), no projection P ^ \ can satisfy 1) 
[P,H] = [P,ha] = and 2) P ha P = Ca P for all a. Then, the non-zero eigenvalues o/ 
have a strictly negative real part. 



Proof: The matrix units \j){k\, j ^ k, are such that 

ho[\j){k\] = -t{E^-E,)\j){k\ . 

Consider hslpjkis)] = ^jk{^) Pjk{^) and the expansion of eigen-matrices Pjki^) = pfk + ^pfk 
and eigenvalues Xjki^) = ~ Ek) +£Vjk to first order in e. Inserted into the eigenvalue 

equation, this yields 



ILo[p5?] +e[u[p'S] + Lo[p«]) ^ -i{E, - Ek)p^ + ^(r],.P$? - ^{E, - Ek)p 



■jk I ' 



whence 

-^[H,pf^] = -^{E, - Ek)p^l , -z[i^ , pja + L,[pS)] = ^,kP^l - i{E, - Ek)pf;i 
Setting p^°^ = one gets 

r,,k = {j\U[\j){mk) = {(jK\j){k\hUk) -l{{j\hiha\j) + {k\hiha\k) 

a 

and, when either i ^ j or k ^ i, 



i{Ei -Ei- Ej + Ek) ■ 
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Therefore, 

a 

< -IY.\^jK\j) - {k\hUk)\' . 

a 

The above inequahty is strict unless {j\hl^ha\j) = \{j\ha\j)\'^ and same for {k\hl^ha\k). Then, 
3fJe(?7jfc) = would imply that, for all a, 

whence, contrary to the assumptions, ha reduces to a scalar multiple of P = + \k){k\ 

on the subspace projected out by P. 

Concerning the perturbation of the eigenvalue 0, choose po such that Lo[p] = 0, but 
D[p] 7^ 0; then, to first order in e, 

{Lo + elD)[p + epi + o{e)] = (^eai + o{e)^ (^p + epi + o{e)^ ^Lo[pi] + D[p] = aip . 

By splitting D[p] = Go o D[p] + Fq o D[p], where Fq = id — Go, one gets the solutions 
Pi = — Lq ^ o Fq o D[p] and D[p] = ap. Since D generates a trace-preserving positive map 
on the 7^^°^ -invariant states, the eigenvalue a must be negative. 

Example 11. Consider D as in ExamplelE, where ha = \ip){a\. Then, [P , \il!){a\] = for 
all 1 < a < d yields P = 1. 



3.1 Entanglement production 

In this section we shall study whether an appropriate, purely dissipative Lindblad dynamics 
can create entanglement even when it is a weak perturbation of a non-entangling unitary 
dynamics. The fact that a Lindblad dynamics that does not include a unitary part is able to 
create entanglement is shown in |1]J in a very concrete example, and the fact that a unitary 
evolution can be added if the invariant state is an eigenstate of the unitary evolution is 
the result in [15]. Here we concentrate on the assumption that this is exactly not the case; 
instead, we tackle the situation where the states invariant under the unitary time evolution 
are all separable. First, we observe that 

Lemma 6. Let the generator Lo be given by a Hamiltonian of the form H = Hi®l + 1®H2 
where Hi has eigenvalues Ei^k o-nd H2 eigenvalues where Ei^k 7^ -£^2,; V/, A;. Then, all 
'-fl^^ -invariant states are separable. If the solutions o/D[p] = 0, where D is as in / f^) . are 
not on the border of the set of separable states, then there exists Eq such that for all e < sq 
the invariant state is unique and separable. 
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Proof: Let sq be smaller than the radius in which the perturbation expansion of p{e) such 
that Le[p(£:)] = converges. Because of Proposition [21 if all solutions of D[p] = are 
invertible, then there can exist only one; as a consequence (see the proof of Lemma [3]), 
there can be only one p{e) within its convergence radius. Further, in a sufficiently small 
neighborhood of a state not on the border of the convex set of separable states, all states 
are separable. 

The following result will instead provide instances of a contrary behavior, more along the 
lines of Example [TJ showing the possibility of creating entanglement by weak dissipative 
perturbations. 

Proposition 3. Consider the generator [2^) with a non-entangling Hamiltonian as in Lemma 
and a dissipative perturbation D such that D[p] =0 has only one solution. Then, there is 
a unique state in the kernel ofh^,, given by the perturbation expansion p{e) = Ylm^^Pn where 

Pn = (-)"((id - o Go o d) o Lo' o d)"[po] . (28) 

Proof: Given a zero-th order approximation po such that Lo[po] = and D[po] = 0, we put 
ourselves in the most general situation where Go o D o Lq ""^ o D[po] 7^ 0. We then add to 
pi = — ILq ^ o D[po] a matrix ai such that Lo[cri] = 0; the new matrix pi = pi + 0"i still solves 
lLo[pi] = — D[po]. Since we want to solve Lo[p2] = — D[pi] by inverting Lq, we seek ai such 
that (see ([24])) 

Go o D[pi] + D[(7i] = . (29) 

Since we assumed that D to have only one state in its kernel, it cannot vanish on Go o D[pi]; 
the latter is a traceless matrix and both its normalized positive and negative parts would be 
states in the kernel of D (see the proof of Lemma [3[). Therefore, 

pi = (^id - D"^ o Go o [pi] = - (^id - D^^ o Go o o Lq ^ o D[po] . 

Iterating this construction, one obtains the contributions to the perturbation expansion as 
in ([28]). 

Remark 2. In general, according to Lemma\^ in order to solve equation [29\) for a generic 
Li in the place of D, one has to consider the map Gi that remains associated to it by the 
time-average and check whether Gi o Go o ]D[pi] = 0. 

Example 12. The map D in Example {TU has only one invariant state; according to ^ its 
inverse is given by D^^[X] = —X on X such that Tr(X) = 0. Therefore, by means of [29\] 
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and [23^] . equation [23^) with n = l and po = ^2^=1 Pa yields 

Pi = -(id - Go o d) o o D[po] = Yl 1^(^)1' - ^^E '^i'^y^^^ (30) 

i=l jjtk 

Pie) = (l+^)El^«l'l^)(^|-^^E^§^^?^l-^')(^l+^(^)^ (31) 

i=l j^k 

Consider the bipartite setting of Lemma\^ and set l<a,P<a, a'^ = d, 

\j) = \a(3) = \a) ® 1/3) where H\a(3) = Eap\a(3) , E^p = + -E2,/3 • 

By transposing the first party with respect to the orthonormal basis {\c()}a=i> in Example^ 
one obtains 

p\e) = {l + e)j^\i;^,\'\a){a\®\P){P\-te J] /^^^ |7)(«I ® + o(^) ■ 

a,/3=l (a,/3)7^(7,5) "'^ 

Suppose ipoiiPi = "^^02/32 = /^^ 7^ '^2 (i^^^ /3i 7^ /32/ i^/ien, choosing an entangled state |0) 
supported only in the subspace spanned by and \a2i32) , one calculates 
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0ai^lV'a2/3i0a2/32V'ai/32 



(0|p^(e)|0)=£ ^ ^ + o{e) . 

r/izs expectation can always be made negative and thus, by applying the partial transposition 
criterion, p{e) results entangled to order e. Notice that the assumptions on the coefficients 
ensure that the projection ofl^p) onto the subspace spanned by |a2/9i) and \aif32) is entangled; 
furthermore, ExamplelTll ensures that p{e) is an asymptotic state for the given time-evolution. 

4 Summary 

A practical control of the invariant states of quantum dynamical semigroups generated by 
Lindblad type generators is still partial, while this knowledge would very much be needed 
due to the fact that quantum information protocols might use suitably engineered open 
quantum time-evolutions to achieve entanglement asymptotically. 

In this paper, we have considered an asymptotic approach to this problem and studied 
the fate of the separable invariant states of two non- interacting quantum systems when their 
Hamiltonian, unitary time-evolution is weakly perturbed by a Lindblad type dissipative 
contribution eLi to the generator Lq. 

The investigation has been conducted by developing the first theoretical steps of a pertur- 
bative approach to the stationary states p{e) of the perturbed generator = Lq -|- e Li and 
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the preliminary results have been applied to show how small suitable dissipative corrections 
to the unitary dynamics, not entangling by itself, may indeed lead to entangled asymptotic 
invariant states, Le[p(£)] = 0. 

The main theoretical tool used to obtain these results has been the practical handling of 
the first steps of an iterative procedure that provides the contributions to the perturbative 
expansion of p{e); the examples we have presented show that in some cases the iterative 
procedure can be performed. From a theoretical point of view, it remains to be understood 
whether the construction of the perturbative contributions to p{e) can always be achieved by 
an appropriate choice of the zeroth order seed as in Section 12.11 or whether the summation 
breaks at a finite order because of lack of analyticity. 
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